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Abstract. The motion of the Earth around the Sun causes an annual change in the mag- 
nitude and direction of the arrival velocity of dark matter particles on Earth, in a way 
analogous to aberration of stellar light. In directional detectors, aberration of weakly in- 
teracting massive particles (WIMPs) modulates the pattern of nuclear recoil directions in 
a way that depends on the orbital velocity of the Earth and the local galactic distribution 
of WIMP velocities. Knowing the former, WIMP aberration can give information on the 
latter, besides being a curious way of confirming the revolution of the Earth and the ex- 
traterrestrial provenance of WIMPs. While observing the full aberration pattern requires 
extremely large exposures, we claim that the annual variation of the mean recoil direction 
or of the event counts over specific solid angles may be detectable with moderately large 
exposures. For example, integrated counts over Galactic hemispheres separated by planes 
perpendicular to Earth's orbit would modulate annually, resulting in Galactic Hemisphere 
Annual Modulations (GHAM) with amplitudes larger than the usual non-directional annual 
modulation. 
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1 Introduction 

Direct dark matter experiments search for energy deposited by the scattering of Weakly 
Interacting Massive Particles (WIMPs) in the dark halo of our galaxy. In non-directional 
direct detection the only dark matter signature is the annual modulation of the interaction 
rate due to the motion of the Earth around the Sun [1] . This motion periodically changes the 
velocity Vj^b of the detector with respect to the Galaxy, or equivalently the average velocity 
of the dark matter particles with respect to the detector. If particles in the dark halo of our 
galaxy are on average at rest with respect to the Galaxy, the average velocity of WIMPs with 
respect to the detector is just — Viab- This is the case of the Standard Halo Model. As in 
this model, we adopt an isotropic Maxwell-Boltzmann (1MB) velocity distribution, but with 
a dispersion and average velocity chosen to reflect our best knowledge of the dark halo of our 
galaxy 

In directional detectors [2], the two dark matter signatures easiest to detect are first 
a departure from isotropy of the recoil directions with respect to the Galaxy, and second 
an average recoil direction in the direction of — Viab [3-8]. Directional detectors, such as 
DRIFT [9], DMTPC [10], NEWAGE [11] and MIMAC [12] use CS2, CF4, ^He or other 
compounds. Refs. [7, 8] find that for an S detector 9 events above 20 keV would be enough to 
reject isotropy, while between 27 and 32 events, depending on the velocity distribution, would 
be enough to confirm the direction of solar motion as the median inverse recoil direction at the 
95% CL. Ref. [7] assumes the recoil directions, including their senses, can be reconstructed 
perfectly in 3d and the background is zero. Also 25 events in the 5 keV to 50 keV range 
in this case were found in Ref. [13] to be necessary to determine the mean recoil direction 
(with a 20° uncertainty), assuming a less optimistic model for a CF4 detector with a 50% 
background contamination. With lOO's of events we could detect a ring of maximum recoil 
rate around — Viab, which would be a secondary indicator of dark matter [14]. 
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After having obtained a dark matter signature, the next goal will be to study the WIMP 
properties, in particular its local velocity distribution. With 100 events Ref. [15] find that 
Viah, the magnitude of Vi^b, and the WIMP velocity dispersion could be determined for a 50 
GeV/c^ WIMP with an uncertainty of about 100 km/s (in CF4). This uncertainty is large 
with respect to the range of the best current estimates of V^ab (between 180 km/s and 340 
km/s; see Table 1) and the dispersion (usually taken to be between 170 km/s and close to 
300 km/s), but such a determination would be of extreme importance. Here we introduce 
the idea of using aberration features of the directional recoil rate to obtain information on 
the local WIMP velocity distribution. 

In directional detectors there will be a change in the recoil direction pattern, which we 
call aberration after a similar effect on the position of stars in the sky. These aberration 
features in a dark matter signal are due to the change in the direction of the arrival velocity 
of WIMPs on Earth caused by Earth's motion around the Sun. These features depend on 
the orbital characteristics of Earth's revolution and on the WIMP velocity distribution, and 
knowing the former they can be used to determine the latter. 

Due to Earth's revolution, the mean recoil direction will change with a period of one 
year. The maximal angular annual separation 7^ between the mean recoil directions in a year 
is inversely proportional to Vi^h- Depending on the velocity of the rotation of the Galaxy 
at the position of the Sun, VcaiRot, this angular separation could be between 18° and 11° 
(for VGaiRot between 180 km/s and 312 km/s in the 1MB model). The percentage error 
in the determination of the WIMP velocity is in this method the percentage error in the 
determination of 7^. Thus we need a precision A-fg of a few degrees. Using the number 
of events given in Refs. [15] and [16] to determine the mean recoil direction within a few 
degrees, in Section 4 we estimate that a number of events between 1400 and a few thousands 
(between 60 kg-yr and 120 kg-yr exposure of CF4) would be needed to measure the angular 
difference in recoil directions integrated over six or three months, which varies in the 1MB 
model between 11° and 7° or 16° and 9°, respectively. 

With more than a thousand events, e.g. more than 60 kg-yr of CF4 (see Section 6), it 
would also be possible to detect the annual modulation of the integrated rate due to aberra- 
tion. The usual non-directional annual modulation of the rate integrated over all directions 
depends only on the annual change in the magnitude of Vjab- The modulation due to aberra- 
tion is much larger. In Section 5 we show that the annual modulation of the rate integrated 
over some Galactic hemispheres, which we call Galactic Hemisphere Annual Modulations 
(GHAM), have much larger amplitudes than the usual non-directional modulation. The 
hemispheres with the largest GHAM are those divided by planes perpendicular to Earth's 
orbit around the Sun, since for them the full orbital velocity is in the direction of one hemi- 
sphere at one time and away from it six months later. The GHAM amplitudes as functions 
of energy depend strongly on the WIMP average velocity and velocity dispersion, and for 
the energy differential rate change sign at much lower energies than the usual non-directional 
annual modulation amplitude (see Section 5). 

We start in Sections 2 and 3 by presenting WIMP and recoil aberration maps. 

2 Annual variation of the WIMP directional flux 

In the 1MB, WIMPs are on average at rest with respect to the Galaxy, and have a Maxwellian 
velocity distribution with dispersion (T„, truncated at the escape speed fosc (with respect to 
the Galaxy). Normalized to 1, the WIMP velocity distribution in the laboratory rest-frame 
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Figure 1. Vsoiar, VEarthRov, and Vlab shown in the Xg-Zg Galactic plane for (a) December 1 and (b) 
May 30. The uncertainties in Vsoiar and VEarthRev are small enough that this result docs not change 
(see text). Here VoaiRot = 180 km/s. 



is given by [17] 
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erf 



vr a 



■ exp 



V 



,,2 -\ 



2al 



(2.2) 



The laboratory is moving with velocity Viab with respect to the Galaxy (thus — Viab is the 
average velocity of the WIMPs with respect to the detector). Following Ref. [18], which gives 
100 km/s and 130 km/s as extreme estimates for the Id velocity dispersion a^jypi^ we take 
o"t, either 173 km/s or 225 km/s. A recent study by the RAVE survey using high velocity 
stars finds an escape speed in the range 498 km/s < Ve&c < 608 km/s at 90% confidence, with 
a median-likelihood value of 544 km/s [19]. We use the median RAVE value, Vesc = 544 km/s 
in this paper. 

The velocity of the lab with respect to the center of the Galaxy Vjab can be divided 
into four components: the Galactic rotation velocity VcaiRot at the position of the Sun or 
Local Standard of Rest (LSR) velocity. Sun's peculiar velocity Vsoiar in the LSR, Earth's 
translational velocity VEarthRev with respect to the Sun, and the velocity of Earth's rotation 
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Date 


VcalRot (km/s) 


IVcalRot + Vsolarl (km/s) 


Viab (km/s) 


December 1 


180 


193 


179.8 


May 30 


180 


193 


208.8 


December 5 


312 


324 


310.5 


June 2 


312 


324 


340.2 



Table 1. Maximum and minimum values of Viab in km/s, for our two choices of the Galactic rotation 
speed VbaiRot • 



around its axis VEarthRot, 

Viab = VcalRot + Vsolar + VEarthRcv + VEarthRot- (2.3) 

The standard value of VcaiRot for the Standard Halo Model is 220 km/s [20]. As dis- 
cussed in Ref. [7], recent studies have found other values for VcaiRot- One analysis found 
VbaiRot = (254 ± 16) km/s [21]. Another study [22] found FcaiRot = (236 ± 11) km/s assum- 
ing a flat rotation curve, while Ref. [23] found values ranging from VcaiRot = (200 ±20) km/s 
to VcaiRot = (279 ± 33) km/s depending on the model used for the rotation curve. We take 
VcaiRot = 180 km/s and 312 km/s as low and high estimates. 

Due to the ellipticity of the Earth's orbit the times of maximum and minimum Vjab 
are not exactly half a year apart. These times depend on VcaiRot- For ^GaiRot = 180 km/s, 
Viab is maximum on May 30 and minimum on December 1. For VcaiRot = 312 km/s, the 
maximum and minimum happen on June 2 and December 5, respectively. Table 1 gives the 
maximum and minimum values of VJab (in km/s) for our two choices of VcalRoti including the 
contribution of the Solar motion Vgoiar- 

Ref. [24] find Vgoiar = (ILll^yg, 12.24^°;^^, 7.25l[]J^) km/s in the Galactic reference 
frame, with additional systematic uncertainties ~ (1,2,0.5) km/s. The uncertainty in Vsoiar 
may shift the dates of the maximum and minimum of V|ab by only one day. Assuming 
VcaiRot = 180 km/s, and taking the low value of Vsoian the maximum and minimum of Vjab 
occur on May 30 and December 2, respectively. For the high value of Vgoiar j the maximum 
and minimum of Vjab occur on May 29 and December 1, respectively. The changed caused 
in the dates of the maximum and minimum Vjab by the uncertainty in V^arthRev is less than 
one day. 

In Fig. 1 we show the components of Vgoiar, VEarthRev, and Viab (VEarthRot is too small 
to be shown) on December 1 and May 30 in the plane perpendicular to VcaiRot (VcaiRot lies 
on the Galactic equatorial plane). This is the Xg-Zg Galactic plane, where the x^-axis points 
towards the Galactic Center and the z^-axis points towards the Galactic North. Notice that 
the component of Viab perpendicular to the Galactic equatorial plane points slightly towards 
the North Galactic Hemisphere (NGH) in December and slightly towards the South Galactic 
Hemisphere (SGH) in May. This is relevant to understand Figs. 2 and 3. 

The plane of the ecliptic (on which VEarthRev lies) is at about 60° of the Galactic 
equatorial plane, so VEarthRev points towards the NGH for half of the year and towards the 
SGH for the other half and VEarthRev — 30 km/s (the average orbital speed is 29.8 km/s) 
is larger than Vgoiar — 18 km/s. The component of Vgoiar in Xg-Zg has a magnitude of 
approximately 13 km/s and a constant direction, while the component of VEarthRev on the 
same plane, has a larger magnitude on May 30 and December 1, about 25 km/s, and inverts 
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Figure 2. Mollweide equal-area projection maps of tlie celestial sphere in Galactic coordinates 
showing the number fraction FwiMp(v,t'q) of m = 100 GeV/c^ WIMPs crossing the Earth per unit 
solid angle as a function of the WIMP velocity direction v with speed larger than Vq = 113 km/s, as 
necessary to produce E}i = 5 keV sulfur recoils. We assume the 1MB with v^sc = 544 km/s, = 173 
km/s and VcaiRot = 180 km/s. (a) (Left panel) on December 1. (b) (Right panel) on May 30. Notice 
the direction of — Viab marked with a cross, (c) Center panel: May— Dec difference of the number 
fractions. The color scale/grayscale shown in the vertical bars corresponds to equal steps between the 
minimum and maximum values in units of sr~^. The directional differential recoil rate at 5 keV in 
CS2 corresponding to this WIMP flux is shown in Fig. 3. 




Figure 3. Directional differential recoil rate at En = 5 kcV in CS2 plotted in Mollweide equal 
area projection maps of the celestial sphere in Galactic coordinates for a WIMP of mass m = 100 
GeV/c^, in the 1MB with VcaiRot = 180 km/s, = 173 km/s and and v^sc = 544 km/s (same 
as in Fig. 2) (a) on December 1, and (b) on May 30. The direction of — Viab is indicated with a 
cross (in the SGH in December and in the NGH in May). The color scale/grayscale shown in the 
vertical bars corresponds to equal steps between the minimum and maximum values given in units of 
10"^ X (po.30-44/kg-day-keV-sr). 



its direction. The result is that the component of Viab in the Xg-Zg plane points North in 
December and is larger in magnitude than in May, when it points South. 

The uncertainties in the measurement of Vsoiar and VEarthRev are small and do not 
change this result. Due to the uncertainty in the ecliptic longitude of the Earth's orbit 
minor axis, Aq = 13° it 1°, there is an uncertainty in T^EarthRev, VEarthRev = (—9.027 it 
0.001, -15.151 ± 0.002, 24.529 ± 0.004) km/s (in the Galactic reference frame) on December 
1 and VEarthRev = (8.979 ± 0.001, 14.725 ± 0.002, -23.798 ± 0.004) km/s on May 30 [25-27]. 
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Figure 4. Difference of tiie directional differential recoil rates at Ef{ = 5 keV in CS2 for a WIMP of 
mass TO = 100 GeV/c^ and (a) May 30-Dec 1, VbaiRot = 180 km/s, cr„ 173 km/s (top left); (b) 
June 2— Dec 5, VcaiRot = 312 km/s, ay = 173 km/s (top right); (c) May 30— Dec 1, VcaiRot = 180 
km/s. (Tj, = 225 km/s (bottom left); (d) June 2— Dec 5, VcaiRot = 312 km/s, (Jy — 225 km/s (bottom 
right). The color scale/grayscale shown in the vertical bars corresponds to equal steps between the 
minimum and maximum values given in units of 10~^ x (po.3(T44/kg-day-keV-sr). 



l^GalRot (km/s) 


180 


312 


Mab (km/s) 


195 


325 


Maximum 


17.8° 


10.6° 


3 months 


16.2° 


9.0° 


6 months 


11.4° 


6.8° 



Table 2. Values of V|ab (in km/s) in March and September, maximal angular separation (March— Sep) 
in degrees, and the angular separation between 3 month averages and 6 month averages. 



The projection of — Viab onto the plane perpendicular to — (VcaiRot + Vsoiar) is maxi- 
mum in March and in September (for VcaiRot > 100 km/s, otherwise they happen at other 
times). In these two months the projection points in opposite directions and Viab happens to 
be the same (see Table 2). Thus the maximal angular separation, 7^ between the directions 
of —Viab in a year is given by the difference between the March and September directions 
(for VbalRot > 100 km/s). 

Table 2 gives the values of VJab (in km/s) in March and September for our two choices 
of VcaiRot) as well as the maximal angular separation 7^ (between the March and Septem- 
ber directions of —Viab), arid the angular separation between the mean recoil directions 
obtained in two three month periods (February to April and August to October) and two six 
month periods (January to June and July to December) centered in March and September, 
respectively. 

The plots of directional rate in the celestial sphere are better viewed in a planar projec- 
tion. Here we use Mollweide equal-area projection maps of the celestial sphere in Galactic 
coordinates. The relationship between the {x,y) coordinates on a Mollweide map and the 
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Galactic longitude and latitude {l,b) is given by (see [28] or [14]) 



2^/2 cos 61 



— TTX 



b = arcsin 



29 + sm{2e) 



vr 



) 



(2.4) 



where 



9 = arcsin 



iy/V2). 



(2.5) 



In Fig. 2 we show a Mollweide projection of the number fraction of 100 GeV/c^ WIMPs 
crossing the detector per unit solid angle as a function of the WIMP velocity direction v 
at the detector. What we plot is the number fraction of WIMPs i^wiMp(v, I'g) moving in 
the direction v with speed higher than Vg. These are WIMPs which can produce a recoil 
momentum of magnitude q or recoil energy Eji = (f' /2M (5 keV in S in Fig. 2) when 
scattering off a nucleus of mass M [29] 



The upper limit of this integral is t;max(v) = -v • Viab + y (v • Viab)^ - Vf^^ + ^2^^, and 
the analytic expression of i^wiMp(v,W5) is given in Eq. 13 of Ref. [29]. The maximum of 
-PwiMp(v, "Wg) happens when v • Vi^b = — ^lab) i-e. in the direction of the average WIMP 
velocity — Viab, i-e. most WIMPs move in the direction opposite to the laboratory motion, 
marked by a cross in the figures. 

The difference between the December and May maps in the upper panels of Fig. 2, given 
in the lower panel of the figure, shows a characteristic aberration pattern with more positive 
numbers in the NGH and more negative ones in the SGH. The number of WIMPs moving 
towards the NGH is larger in May than in December, while the opposite is true for WIMPs 
moving towards the SGH. This is what we would expect from the change in the direction of 
—Viab from the NGH in May to the SGH in December. 

3 The directional differential recoil spectrum 

The aberration features in directional detection depend only on the Radon transform of the 
WIMP velocity distribution, and not on the particular type of WIMP-nucleus interaction. 
The directional differential recoil spectrum as a function of the recoil momentum q is given 
in terms of the 3-dimensional Radon transform of the WIMP velocity distribution /lab [5] as 



Here p is the dark matter density in the solar neighborhood and m is the WIMP mass. 
The sum is over the nuclear species i in the target, and Cj and /ij are the mass fraction 
and the reduced WIMP-nucleus mass for nuclide z, respectively, d^q = d(/)dcos9 denotes an 
infinitesimal solid angle around the recoil direction q = q/q, q = |q| is the magnitude of the 
recoil momentum, g/2/i = Vq is the minimum velocity a WIMP must have to impart a recoil 
momentum q to the nucleus of mass M, or equivalently to deposit a recoil energy Eji = 
q^ /2M, /i = mM/{m + M) is the reduced WIMP-nucleus mass, iTj(g) is the WIMP-nucleus 
scattering cross section which can be split into spin-independent (SI) and spin-dependent 




(2.6) 





(3.1) 



(SD) parts, cri{q) 



af\q)+af^{q). 
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The Radon transform in the laboratory frame for the truncated Maxwelhan WIMP 
velocity distribution in Eq. 2.1 is [17] 



lab 



2/1 



iVesc(2vra2)l/2 



exp 



[(g/2/i) + q • V: 



labj 



2al 



exp 



2al 



(3.2) 



if {q/2jj,) + q • Vjab < Vcsc, and zero otherwise. 

The recoil momentum q is measured in a reference frame fixed to the detector. The 
detector frame is at some orientation in the laboratory frame, which we define as fixed to 
the Earth with axes pointing to the North, the West and the Zenith. The transformation 
equations between the detector frame and the laboratory frame can be conveniently written in 
terms of direction cosines measurable in any experiment. This is formulated in the Appendix 
A of Ref. [14], where we also give the transformations from the laboratory frame in any 
location on Earth to the Galactic reference frame. The transformations in Appendix A of 
Ref. [14] take into account Earth's rotation around its axis, which is usually neglected. We 
use Galactic coordinates, but one could use any coordinate system to describe the aberration 
features. 

To give concrete examples of the GHAM we choose WIMPs with SI interactions and 
equal couplings to protons and neutrons interacting in a CS2 detector or with SD interactions 
in CF4. For these WIMPs, the SI cross section in Eq. 3.1 is crf^iq) = iJ.fAfapSi{q)/iJ,p, where 
Hp = mnip/ {m + nip) is the WIMP-proton reduced mass, Up is the SI WIMP-proton cross 
section, Ai is the mass number of the nuclear species i, and Si{q) is the nuclear form factor, 
for which we use the Helm [30] nuclear form factor normalized to 1. The SI directional 
differential recoil rate is therefore [5] 



1.306 X 10' 



events 



kg-day-keV-sr 



47rm/i^ 



(3.3) 



Here po.3 is the dark matter density in units of 0.3 GeV/c^/cm'^, CT44 is the WIMP-proton 
cross section in units of 10~^^ cm^, /ip and m are in GeV/c^, and /lab is in (km/s)~^. 

The SD directional differential recoil rate in the proton-odd approximation in CF4 is [5] 



dEn dflq 



130.6 



events 



kg-day-keV-sr 



5(0.647)C./,..^^^^„ 



(3.4) 



where cr^^g is the spin-dependent cross section off a proton in units of 10"'^^ cm^ (consistent 
with experimental bounds from direct detection and neutrino experiments [31]), and Cp = 
AMf/[AMf + Mc) is the mass fraction of F in CF4. Note that while the F nucleus has spin 
^, the C nucleus has no spin. 

The standard value of the local dark matter density \s p = 0.3 GeV/c^/cm'^, and this 
is what we use here. However, one should keep in mind the large uncertainties in this 
parameter. Recent astronomical constraints are consistent with 0.2 GeV/c^/cm"^ < p < 0.4 
GeV/c^/cm'^ for a spherical dark matter halo profile, and up to 20% larger for non-spherical 
haloes [32]. Using a halo model independent method, Ref. [33] finds p = 0.43 =b 0.11 it 0.10 
GeV/c^/cm^, with uncertainties from two different sources. Ref. [34] finds the density for 
a specific simulated galaxy resembling the Milky Way is 21% larger than the mean value of 
p = 0.39 GeV/c^/cm^ obtained in a previous study [35] in which spherical symmetry was 
assumed. 
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Fig. 3 shows the directional differential recoil rate (given in Eq. 3.3) in CS2 at = 5 
keV on December 1 (left panel (a)), and May 30 (right panel (b)), the times at which Vjab 
is maximum or minimum, plotted in Mollweide maps of the recoil direction q in Galactic 
coordinates for a WIMP of mass m = 100 GeV/c^. The 1MB values in this fi gure are 
VcaiRot = 180 km/s, £7^ = 173 km/s, and Ucsc = 544 km/s. The direction of the average 
WIMP velocity with respect to the detector, — Viab, is indicated on the maps with a cross. 
The difference between the two Fig. 3 maps is shown in the upper left panel of Fig. 4 (in 
units of 10~^ X (/3o.3<744/kg-day-keV-sr)) together with similar plots obtained by changing 
only Flab (i-e. VcalRot) and to their maximum and minimum values. 

As seen in the four panels of Fig. 4, increasing the value of VcaiRot results in a larger 
anisotropy in the aberration features. This is expected since the anisotropy in the flux of 
WIMPs arriving on Earth and thus in the differential recoil rate increases with VQaiRot (the 
flux would be isotropic if VcaiRot = 0). On the other hand, increasing makes the difference 
maps more isotropic, also as expected (cj^ = would mean that all WIMPs arrive from only 
one direction). The most anisotropic aberration features (for the VcaiRot and values that 
we considered) is shown in the upper right panel of Fig. 4 in which VcaiRot is the largest and 
ay the smallest, 312 km/s and 173 km/s, respectively. Notice that the right panels of Fig. 4 
(with VcaiRot = 312 km/s) show a more negative rate difference (darker region) in the NGH 
compared to the left panels. This is due to the higher visibility of the ring-like feature in the 
recoil maps of June and December for larger VcaiRot (see Figs. 5.b and 6 of Ref. [14] in which 
the ring is more visible for VcaiRot = 312 km/s). Changing the value of t^esc does not change 
the maps significantly, and thus we show the results only for one. 

Examining the four panels of Fig. 4, we can see that the positive directional differential 
rate differences (lighter regions) are mostly in the NGH and the negative rate differences 
(darker regions) are mostly in the SGH. This effect persists in the energy and time integrated 
rates. Using the same parameters of Fig. 3, Fig. 5. a shows the December and May difference 
in the rates averaged over the 5 keV to 20 keV energy range. Fig. 5.b shows the difference 
in the directional differential rates averaged over three months (December, January, and 
February period with respect to the June, July, and August period), and Fig. 5.c shows the 
difference in the time averaged (over the same three months) and energy averaged (in the 5 
keV to 20 keV range) of the same rates. Clearly, the aberration features persist in the three 
panels of Fig. 5. 

In experiments that are not able to distinguish the sense of a nuclear recoil (such as 
NEWAGE), the directional differential recoil rate detected in opposite Galactic quadrants 
are identical. Fig. 6 shows the directional recoil rate without sense discrimination in CS2 
at En = 5 keV plotted in Mollweide maps of the recoil direction q in Galactic coordinates 
for a WIMP of mass m = 100 GeV/c^, and assuming the 1MB with Vesc = 544 km/s, 
^aiRot = 180 km/s and = 173 km/s. The directional differential rate is given on December 
1 in the left panel (a), and on May 30 in the right panel (b). The center panel (c) gives 
the difference in the directional differential rate between May and December, in units of 
10~^ X (/3o.3<T44/kg-day-keV-sr). Having sense discrimination in 3d makes the detection of 
the aberration features easier and gives the largest annual modulation amplitudes. Thus in 
the rest of this paper we give the results assuming the directional detector has 3d capability 
of reconstruction of recoil directions, including their senses. 

Even with sense discrimination, the aberration patterns require a very large number of 
events to be observed. We concentrate in the following sections on aberration features easier 
to detect. 
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Figure 5. Same parameters as in Figs. 3 and 4. a. (a) (Left) (May 30) — (Dec 1) difference of the 
directional differential recoil rates averaged over an energy interval between 5 keV and 20 keV, (b) 
(right) difference of the time averaged directional differential recoil rates over three months at 5 keV, 
and (c) (center) difference of the directional differential recoil rates averaged over three months and 
over the 5 keV to 20 keV energy interval. In (b) and (c) the average rate of December, January, and 
February is subtracted from the average rate of June, July, and August. 




Figure 6. (a) and (b) same as Figs. 3. a and 3.b and (c) same as Fig. 4. a but without sense discrimi- 
nation. 



4 Maximal angular variation of the mean recoil direction 

The largest uncertainty in Vj^b is due to our uncertain knowledge of the velocity of the 
rotation of the Galaxy at the location of the Sun. Thus we can use the change in the average 
velocity of WIMPs with respect to Earth in a year to measure VcaiRot • A simple and rather 
accurate evaluation of the maximal angular separation 7^ between the directions of — Viab in 
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a year is 



7s ^ 



2^EarthRev sill 60° 
^GalRot 



16.5 



O 



( 



180 km/s 



^GalRot 



) 



(4.1) 



which gives from 16.5° to 9.5° for VcaiRot between 180 km/s and 312 km/s. The maximal 
angular separation, 7^ is the difference between the — Vjab directions in March and September. 
This angle 7^ in degrees is plotted in Fig. 7, together with the angular separations between the 
mean recoil directions obtained in two six month periods or two three month periods centered 
in March and September, respectively. 7^ goes between 16.2° and 9.0° for the three months 
averaged directions, and 11.4° and 6.8° for the directions averaged over six months (see 
Table 2). The error AVcaiRot in the determination of VcaiRot using 7^ can be read off of Fig. 7 
(see below). A simple estimate, 7^ ~ ^0^11101 indicates that AVcaiRot/VGaiRot - ^Is/ls- 
Thus we need A7S to be a few degrees at most. 

Refs. [16] and [6] give estimates of the number of events needed to measure the mean 
recoil direction with an error of a few degrees in two different cases. Ref. [16] finds that with 
an exposure of 30 kg-yr in CF4 the mean recoil direction could be determined with an error 
of 2.5°. For a 50 GeV/c^ WIMP with SD interactions and af^ = 10"^^ cm^ and recoils 
with energy between 5 and 50 keV, as assumed in Ref. [16], 30 kg-yr corresponds to 660 
events in this energy range. In Ref. [16], Vjab is fixed to 220 km/s, however it is reasonable 
to assume that with a variable VJab the mean recoil direction would be determined with a 
similar error of a few degrees. If the data would be split in two sets (January to June and 
July to December), the six months averaged directional separation 7^ could be determined 
with double the number of events and exposure, i.e. 1400 and 60 kg-yr. The three months 
averaged 7^ would require four times the total number of events and exposure, i.e. 2800 and 
120 kg-yr. Assuming an error of A7, ~ \/2 X 2.5° ~ 3.5°, it can be seen from Fig. 7 that the 
measur ed VbaiRot would be 180l^[] km/s or 312^^^° km/s using the three month averaged 
(green) curve, and 180j]4Q km/s or 3121;^^q km/s using the six month averaged (orange) 



Ref. [6] shows in its Fig. 1 that more than 10"^ events above 20 keV would be necessary in 
S to determine the mean recoil direction with an error smaller than 5° (in this case m = 100 



Thus, both Refs. [16] and [6] imply that a few thousand events would be necessary to 
measure 7^ . With this number of events it might be possible to detect the annual modulation 
of the rate integrated over Galactic hemispheres to which we devote the following two sections. 

We would like to point out that VcaiRot could also be obtained from the annual average 
value of the projections of recoil momenta onto planes containing the mean recoil direction. 
One such plane could be chosen to be perpendicular to a vector at 90° from both Earth's 
orbital plane and the mean recoil direction. Another could be the plane perpendicular to 
this one and also containing the mean recoil direction. The average of the projections of the 
recoil momenta onto these two planes should be ~ VcaiRot while the average of the projections 
onto the plane perpendicular to the two just mentioned would be zero. We do not develop 
this method any further for the time being, because we think it would not require a smaller 
number of events than the previous method. 

5 Annual modulation integrated over direction 



curve. 



GeV/c2). 



The annual modulation of the differential rate is due to the change in the magnitude of 
T^ab during a year [1]. At high energies, the differential rate is maximum and minimum 




Figure 7. Maximal angular separation 7s between the — Viab directions in a year (black line), and 
between the mean recoil directions obtained in two three month periods (green line) and two six 
month periods (orange line) , plotted as a function of VcaiRot • 



respectively when Vjab is maximum and minimum, but at low energies the phase is inverted 
and the differential rate is maximum when V\ah is minimum and vice versa. 

In order to compute the annual modulation amplitude of the energy differential rate, 
we integrate the directional differential recoil rate over direction, and define 



f I dRraax dR 

VdE^ J ~ dEn dEn ~ J [dEn dn^ ~ dEn dn, ' ^ ' ' 



where the subscripts "max" and "min" refer to the maximum and minimum differential rates 
during a year. In the plots we show the amplitude at high energies as positive (thus it 
becomes negative at low energies). 

In directional detection, we can integrate over the recoil directions pointing to half of 
the sky and define the annual modulation of the rate so obtained. This is what we call 
the Galactic Hemisphere Annual Modulation (GHAM). Choosing well the hemisphere, the 
GHAM amplitude is larger than the usual annual modulation amplitude, thus easier to detect. 

5.1 Maximum GHAM amplitudes 

As mentioned earlier, in the Mollweide figures showing the difference in the directional dif- 
ferential recoil rates between the two dates at which Viab is maximum and minimum, one can 
see that the positive rate differences are mostly in the NGH and the negative rate differences 
are mostly in the SGH. This happens because the plane of the Earth's orbit around the Sun 
is at about 60° of the Galactic equatorial plane. We can compute the annual modulation 
amplitude for any Galactic hemisphere, but those with the largest GHAM are divided by 
planes that are perpendicular to the Earth's orbit around the Sun, since then the total or- 
bital velocity is in the direction of one hemisphere at some time and away from it half a year 
later. 

Any hemisphere can be defined by a vector perpendicular to a plane dividing the sky in 
two, and pointing towards the particular hemisphere. We call this vector the hemisphere's 
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Figure 8. Annual modulation amplitude (Eq. 5.1) of the S energy differential recoil rate integrated 
over different Galactic hemispheres (colored lines) and the whole sky (black line) as a function of 
recoil energy En for m = 100 GeV/c^, Wosc = 544 km/s and (a) VcaiRot = 180 km/s and ~ 173 
km/s (top left), (b) VGaiRot ~ 312 km/s and ~ 173 km/s (top right), (c) VcaiRot = 180 km/s and 
(Tu = 225 km/s (bottom left), and (d) VcaiRot = 312 km/s and a^, ~ 225 km/s (bottom right). The 
red, purple, cyan, and blue curves are for different dates (rate at Dl minus rate at D2) in a year 
separated by three months: Dl January 1, April 1, July 1, and October 1, respectively, and D2 July 
2, October 4, December 31, and March 29, respectively. The green and orange lines show the GHAM 
with largest amplitude at low and high energies, respectively, and they correspond to complimentary 
hemispheres. For the green lines, the Dl— D2 dates are Dec 1— May 30 in (a) and (c), and Dec 5— June 
2 in (b) and (d). For the orange lines the Dl and D2 dates are exchanged. 



pole. Let us consider the Galactic hemisphere with the vector VEarthRov at a particular time 
of the year, the date Dl, as its pole. The maximum differential rate for this hemisphere 
happens on the date Dl, and the minimum rate happens on the date D2 when VEarthRev 
points in the opposite direction. 

The two hemispheres with the maximum GHAM amplitudes at high and at low ener- 
gies have as their poles VEarthRov on the dates at which Vjab is maximum and minimum, 
respectively. 

Figs. 8 and 9 show the annual modulation amplitude t\{dRl dE^) of S recoil rates as 
a function of the recoil energy E'/j for m = 100 GeV/c^ and m = 7 GeV/c'^, respectively 
integrated over different hemispheres (colored lines) and over all directions (black line) . In the 
four panels of the figure we use one of the four different combinations of VcaiRot = 180 km/s 
and 312 km/s, and = 173 km/s and 225 km/s. The escape velocity is always Vqsc = 544 
km/s. The red, purple, cyan, and blue curves are for hemispheres having as pole VEarthRev 
on different dates in a year separated by three months. The green and orange curves show the 
GHAM with largest amplitude at low and high energies, respectively, and they correspond 
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Figure 9. Same as Fig. 8 but for m = 7 GeV/c^. 



to complimentary hemispheres. For the green hnes, the Dl— D2 dates are Dec 1— May 30 in 
the (a) and (c) panels, and Dec 5— June 2 in the (b) and (d) panels. For the orange lines 
the Dl and D2 dates are exchanged. The annual modulation amplitude (black line) is thus 
the orange line minus the green line, and goes to zero when they are equal. Notice that the 
GHAM amplitudes are much larger and become negative at much smaller energies than the 
usual non-directional annual modulation amplitude. 

In Figs. 8 and 9, one can see that the annual modulation amplitudes strongly depend on 
V|ab and (Tt,, although more work is necessary to quantify this dependence. For smaller V|ab, 
the annual modulation amplitude, which is approximately (VEarthRev/Mab) times the annual 
average rate, is larger. A smaller also leads to a more anisotropic recoil rate, which results 
in a larger rate difference. Thus the largest modulation amplitudes happen for the smallest 

and a"„. In Figs. 8 and 9, the maximum amplitude in the (a) panels are about a factor 
of two larger than the maximum amplitudes for the largest Vjab and a^,, shown in the (d) 
panels. Not only the magnitude of the amplitudes but their shape as a function of En change 
with V\sh and Uy. As the energy increases (and thus Vq increases), for larger Vjab or o"^, there 
are more WIMPs with velocity larger than Vq, and thus there is a larger fraction of recoil 
events at higher energies. Therefore, increasing V^ab or f^v^ flattens the recoil spectrum and 
thus also the rate difference. This can be seen in Figs. 8 and 9 where in panels (a) the slope 
of the curves beyond the peak is the largest. 

Measuring the energy differential rate requires very large statistics. Thus we explore 
next the GHAM in the energy-integrated rate. Let i?max = f^^ dEji{dR^3_x/dEii) and 

-Rmin = dEji{dRraiii/ dEj{) denote the maximum and minimum energy-integrated rates 
during a year integrated over a particular energy interval between Ei and E2. The annual 
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Figure 10. Same as Fig. 8 but sliowing tlie annual modulation amplitude of the energy-integrated 
recoil rate as a function of the recoil energy Efj and integrated above Efi and over the two hemispheres 
with the GHAM with largest amplitude at high (orange curve) and low (green curve) energies, as well 
as over the total sky (black curve). 



modulation amplitude of the energy-integrated rate is then 

Ai? = iJ^ax - i?min = A i^-^j cIEr. (5.2) 

Fig. 10 shows the annual modulation amplitudes of S recoil rates for m = 100 GeV/c^ 
integrated above Er, as a function of Er, for the two hemispheres with the GHAM with 
largest amplitude at high (orange curve) and low (green curve) energies, as well as over 
the total sky (black curve). As explained before, these two are complimentary hemispheres: 
in one of them the rate is maximum when V|ab is maximum and minimum when V^ab is 
minimum (orange line), and in the other the opposite is true (green line). In Fig. 10, we 
show the difference between the rates at date Dl and date D2, choosing Dl as the date at 
which Viab is maximum and D2 as the date at which Viab is minimum. Thus the orange 
GHAM amplitude is positive, the green is negative, and the usual non-directional annual 
modulation is the sum of the two. The orange GHAM amplitude reaches values between two 
and three times larger than the usual non-directional modulation amplitude, at low energies. 
The parameters used in the four panels are the same as the parameters used in Fig. 8. 

In Fig. 10 the GHAM with largest amplitude at high energies (orange curve) is larger 
in magnitude than that at low energies (green curve). This is because the area under the 
orange curve in Fig. 8 (which shows the modulation amplitude of the energy differential recoil 
rate) is larger than the area under the green curve in the same figure. The dependence of 
the energy-integrated annual modulation amplitudes on the values of V\ah and ay is similar 
to that of the energy differential rate (Fig. 8). The peak of the largest GHAM amplitude is 
largest for the smallest Viab and (T„. 
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Figure 11. Annual modulation amplitude (NGHAM, SGHAM and total) of the S energy differential 



recoil rate, Eq. 5.1, as a function of the recoil energy En for m = 100 GeV/c 



544 km/s and 



(a) VcaiRot = 180 km/s and ct^ = 173 km/s (left), and (b) VcaiRot = 312 km/s and ~ 225 km/s 
(right). The black, green and orange curves arc for the total (i.e. non-directional), NGHAM, and 
SGHAM amplitudes, respectively. Notice that the GHAM amplitudes are much larger than the usual 
non-directional annual modulation amplitude and change sign at lower energies. 



5.2 NGHAM and SGHAM amplitudes 

In hemispheres other than those defined in Section 5.1, the GHAM are smaller. As an ex- 
ample, we consider here the North and South Galactic hemispheres, i.e. the North GHAM 
(NGHAM) and the South GHAM (SGHAM). The amplitudes of these modulations are re- 
spectively maximum when the component of VEarthRov towards either the North or South 
directions (i.e. in the positive and negative Zg direction) is maximum, which happens on 
December 22 and June 21, respectively. These amplitudes and the total (i.e. for the whole 
sky) modulation amplitudes for the S energy differential rate (Eq. 5.1) are shown is Fig. 11 
as a function of the recoil energy Er, for m = 100 GeV/c^, Vesc = 544 km/s, and extreme 
combinations of VcaiRot and a^, 180 km/s and 173 km/s in the left panel, and 312 km/s and 
225 km/s in the right. As explained earlier, the combination of small VcjaiRot and o'v results 
in the higher modulation amplitudes, whereas the opposite happens for the combination of 
large VcaiRot and a^. 

We find that both GHAM amplitudes in Fig. 11 are smaller than the GHAM amplitudes 
of the previous section (Fig. 8), but still larger and change sign at much lower energies than 
the total sky amplitude. The NGHAM amplitude is positive except at very low energies and 
the SGHAM amplitude is negative, having chosen to present the difference in rates at the 
same two dates. This is what one expects because the mean recoil direction, the direction of 
~Viab, points slightly North in June and slightly South in December (see Figs. 1, 2 and 3). 
Thus going from June to December, the recoil rate decreases in the NGH and increases in 
the SGH. 

6 Statistical test of detectability of the GHAM 

Here we estimate the minimum number of events necessary to detect the GHAM. Let 
and Ni denote the number of recoil events during the annual half-cycles with high, and low 
rates respectively, in a particular recoil energy interval between Ei and £'2- The number of 
events in both cycles is = Nh + Ni. For the annual modulation to be observable at the 3a 
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level we require 



Nh-Ni>3a^ 30V72, (6.1) 



where a ~ N/2 is the standard deviation of or Ni, given that Nfi ~ Ni for a small 
annual modulation. The energy-integrated rate can be written as R{t) = Rq + Rm cos{ujt) 
with uj = 27r/(l yr). The unmodulated part is Rq = {Rma.x + -Rmin)/2 and the modulated 
part has amplitude Rm = {Rma.x — -Rmin)/2. If the detector exposure is MT, then Nh = 
MT{Rq/2 + i?m/7r) and Ni = MT{Rq/2 - Rm/-^)- The factors of l/vr come from integrating 
the modulated part of the rate, Rm cos{ujt) over a half-cycle. 

We define the relative modulation amplitude A (taking into account the signal only and 
neglecting background) in terms of the maximum and minimum recoil rates integrated over 
energy and direction, 

A Rmax -^min Rm /„ „\ 

Tr ^~R~- ^^■^> 

-n-max T Ji-min -'''0 

Thus, Eq. 6.1 becomes 

iV > ^ = iVmin, (6.3) 

where k = vr^/S for a sinusoidal modulation (it would be slightly different otherwise). Alter- 
natively, Eq. 6.3 can be written as 

9k 18k 
A'^Ro ~ A^{R 

max ~l~ Rmin) 

We can apply this condition to any annual modulation. The condition for the ob- 
servability of the GHAM is Nh > 9k/ Aj^, where Nh is the number of events towards one 
particular Galactic hemisphere H, and Ah is the relative modulation amplitude for that 
Galactic hemisphere. One can also find the particular exposure {MT)H,min = 9'^/(^//-Ro,h) 
for the hemisphere H, where Ro,h is the average rate in the hemisphere H. Using the average 
rate Rq^h' in the complimentary hemisphere H', one can find the minimum number of events 
needed with {MT)H,min in the complimentary hemisphere H', Nh'. min = ^H,ramRo,H' / Ro,H ■ 
Thus the minimum number of events observed in all directions necessary to detect the GHAM 
in the initial hemisphere H is 

AtGHAM ^ ^^^^.^ ^ NH',min. (6.5) 

To see if the GHAM would be easier to detect compared to the total annual modulation, 
we need to compare the minimum number of events A^min necessary to detect the total annual 
modulation and A^?^^ in all directions necessary to detect the GHAM. 

Fig. 12 shows the minimum number of S recoils detected (in all directions) with energy 
larger than Eji, Eqs. 6.3 and 6.5, needed to detect the total annual modulation (black curve), 
and the GHAM with the largest amplitude at low (green curve) and high (orange curve) 
energies. In the four panels of the figure we use one of the four different combinations of 
^aiRot = 180 km/s and 312 km/s, and = 173 km/s and 225 km/s, for m = 100 GeV/c^ 
and Vesc = 544 km/s. Fig. 12 shows that one needs a minimum of a few thousand events to 
observe the largest GHAM amplitudes, and one would need at least 10 times more to observe 
the usual non-directional annual modulation. 

Notice that although Amin is smaller at larger energies, the average rate Rq decreases 
with energy, thus the necessary exposure increases at higher energies. This can be seen in 
Fig. 13 which shows the minimum exposure (Eq. 6.4) as function of Eji needed to detect the 



-17- 



S, m=100 GeV, VGaiRot=180 km/s, o-,,=173 km/s 



S, m=100 GeV, V'GaiRot=312 km/s, o-,,=173 km/s 



1x10^ 
SxlO"* 

§ 1 X 10^ 
A 5000 




1000 
500 



1 X 10, 
5x10"^ 

§ 1x10 
A 5000 



1000 
500 

100 




Er (keV) 

S, m=100 GeV, yGaiRot=180 km/s, o-„=225 km/s 



Er (keV) 

S, m=100 GeV, VGaiRot=312 km/s, o-,,=225 km/s 



1x10- 
5x10' 

^1 1 X 10"* 
A 5000 




1000 
500 



40 60 
Er (keV) 



40 60 
Er (keV) 



Figure 12. Minimum mimbcr of events with energy larger than Ef( (Eqs. 6.3 and 6.5) in ah directions 
needed to detect the usual non-directional annual modulation (black curve) and the GHAM with 
largest amplitude at low (green curve) and high (orange curve) energies in an S detector. The 
parameters used are m = 100 GeV/c^, Vesc = 544 km/s and (a) VcaiRot = 180 km/s and ct„ = 173 
km/s (top left), (b) VcaiRot = 312 km/s and = 173 km/s (top right), (c) VcaiRot — 180 km/s and 
(Tt, = 225 km/s (bottom left), and (d) VcaiRot = 312 km/s and CTu = 225 km/s (bottom right). 



non-directional annual modulation (black curve), and the GHAM with largest amplitude at 
low (green curve) and high (orange curve) energies in the rate integrated above i?/? for an S 
detector. The parameters used are the same as those in Fig. 12. Because increases and 
i?o decreases with energy, the increase or decrease of (MT)jnin with energy depends on the 
rate of change with energy of and Rq. Thus (Mr)miii has a minimum at a certain recoil 
energy. 

Fig. 14. a shows the minimum number of events in all directions, with energy above Er, 
needed to detect the GHAM with largest amplitude at high energies in a CF4 detector for a 
50 GeV/c'^ WIMP with SD interactions. We can use SD cross sections larger than the SI cross 
sections = 10~^^ cm^) without violating experimental bounds, which results in smaller 
exposures needed to detect the GHAM. Fig. 14. b shows the minimum exposure needed to 
detect the GHAM with largest amplitude at high energies in the SD rate integrated above 
in CF4. The four curves in Fig. 14 are for the four different combinations of VcalRot = 180 
km/s and 312 km/s, and = 173 km/s and 225 km/s. Fig. 14. a shows that one needs 
between a few hundred and a few thousand events depending on the value of VcaiRot and cr„ 
and the energy interval, to observe the largest GHAM amplitude. From Fig. 14. b one can see 
that the minimum exposure needed to detect the largest GHAM in the SD rate is between 
0.06 ton-yr and a few ton-yr. Figs. 12 and 14. a show that the combination of low VcaiRot 
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Figure 13. Exposure MT (in ton-yr) corresponding to the minimum number of events of Fig. 12 
needed to detect the usual non-directional annual modulation (black curve) and the GHAM with 
largest amplitude at low (green curve) and high (orange curve) energies, assuming p = 0.3 GeV/c^/cm'^ 
and (Tr, = 10^"*^ cm^. 




Figure 14. (a) (Left) Minimum number of events (in all directions) and (b) (right) minimum exposure 
MT (in ton-yr) needed to detect the GHAM with largest amplitude at high energies in the SD rate 
integrated above in CF4. The different curves in each plot are for different values of VcaiRot and 
(Tt,, 180 km/s and 173 km/s (thick solid), 312 km/s and 173 km/s (dashed), 180 km/s and 225 km/s 
(dotted), and 312 km/s and 225 km/s (thin solid), respectively. In (b) we assume p = 0.3 GcV/c^/cm'^ 
and cr^D = 10-39 cm2. 



and gives the smallest minimum number of events. 



- 19 - 



7 Anisotropic Logarithmic-ellipsoidal models 

So far we have used the isotropic Maxwell-Boltzmann model. While the details of the velocity 
distribution are important for the detectability of aberration effects, their existence depends 
only on the motion of the Earth around the Sun, and thus the aberration effects are general 
for any WIMP distribution. A complete analysis of aberration effects in non-Maxwellian or 
non-isotropic distributions is beyond the scope of this paper. However, in Fig. 15 we present 
one such example: an anisotropic logarithmic-ellipsoidal model of Ref. [36] (discussed in 
Section 5 and Appendix B of Ref. [14]). The equations for the velocity dispersions in the 
direction of the Galactic center, the Galactic rotation, and the north Galactic pole, axg, Cyg, 
and (Tzg, as well as the constants p and q used to describe the axis ratios of the density 
ellipsoid, and the constant F used to describe the velocity anisotropy, are given in Appendix 
B of Ref. [14]. In the example presented in Fig. 15, the Solar system is on the minor axis 
with p = 0.72, q = 0.7, and F = 4.02. For VcaiRot = 180 km/s, this choice of parameters 
corresponds to axg = 134 km/s, ayg = 32 km/s, and azg = 177 km/s, and for VcaiRot = 312 
km/s, it corresponds to axg = 233 km/s, ayg = 56 km/s, and a^g = 307 km/s. In Fig. 15 
m = 100 GeV/c^, and the solid and dashed curves are for VcaiRoti 180 km/s and 312 km/s, 
respectively. 

Fig. 15. a shows the annual modulation amplitude A{dR/dEii) of S recoil rates as a 
function of integrated over the two hemispheres with the GHAM with largest amplitude 
at high (orange curve) and low (green curve) energies, as well as over the total sky (black 
curve). For the solid and dashed green curves, the Dl— D2 dates are Dec 1— May 30 and 
Dec 5— June 2, respectively. For the orange curves the Dl and D2 dates are exchanged. 
Fig. 15. b shows the annual modulation amplitudes integrated above Er, as a function of Er, 
for the two hemispheres mentioned and over the total sky. In this figure, the solid curves 
are for May 30— Dec 1, and the dashed curves are for June 2— Dec 5. Figs. 15. c and 15. d 
show the minimum number of S recoils detected in all directions, and the minimum exposure, 
respectively, needed to detect the total annual modulation (black curves), and the GHAM 
with the largest amplitude at low (green curves) and high (orange curves) energies in the 
rate integrated above Er. 

One can see that the magnitude and shape of the annual modulation amplitudes as a 
function of Er in the anisotropic model assumed (Figs. 15. a and 15. b) are quite different 
from those in the 1MB (Figs. 8 and 10). In this anisotropic logarithmic-ellipsoidal model, the 
velocity dispersions depend on VcaiRoti and thus the GHAM amplitudes strongly depend on 
I^ab- Similar to the 1MB, the largest modulation amplitudes happen for the smallest T^abi 
and increasing Mab) flattens the recoil spectrum and thus also the rate difference. Fig. 15.c 
and 15. d show that one needs a minimum of a few hundred events corresponding to at least 
a few ton-yr of S to observe the largest GHAM amplitudes in this anisotropic model. This 
number of events and the corresponding exposure are of the same order of magnitude as for 
the isotropic Maxwellian case. 

8 Summary 

In this paper we study changes in the recoil direction pattern in directional detectors, which 
we call aberration after a similar effect on the position of stars in the sky. These aberration 
features in a dark matter signal are due to the change in the direction of the arrival velocity 
of WIMPs on Earth, caused by Earth's motion around the Sun. These features depend on 
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Figure 15. An example of an anisotropic logarithmic-ellipsoidal model with p = 0.72, q = 0.7, 
r = 4.02, and m = 100 GeV/c^, VcaiRot = 180 km/s (sohd curves), and VcaiRot = 312 km/s (dashed 
curves), (a) (Top left) Annual modulation amplitude of the S energy differential rate integrated over 
the two hemispheres with the GHAM with largest amplitude at high (orange curves) and low (green 
curves) energies, as well as over the total sky (black curves). For the solid and dashed green curves, 
the Dl— D2 dates are Dec 1— May 30 and Dec 5— June 2, respectively. For the orange curves the 
Dl and D2 dates are exchanged, (b) (Top right) Same as (a) but showing the annual modulation 
amplitude of the energy-integrated rate as a function of the recoil energy En and integrated above 
Eji and over the same two hemispheres, as well as over the total sky. The solid curves are for May 
30— Dec 1, and the dashed curves are for June 2— Dec 5. (c) (Bottom left) Minimum number of events 
(in all directions) and (d) (bottom right) minimum exposure AIT (in ton-yr) needed to detect the 
total annual modulation (black curves), and the GHAM with the largest amplitude at low (green 
curves) and high (orange curves) energies in the rate integrated above En in S. In (d) we assume 
p = 0.3 GeV/cVcm^ and dp = 10" 



cm2 



the orbital characteristics of the Earth's revolution and on the WIMP velocity distribution. 
Knowing the former, aberration features can be used to determine the latter, besides being 
a curious way of observing Earth's revolution and confirming the galactic provenance of 
WIMPs. While observing the full aberration pattern would require extremely large exposures, 
we show that the annual variation of the mean recoil direction or of the event counts over 
specific solid angles may be detectable with moderately large exposures. 

The maximal angular annual separation 7^ between the mean recoil directions in a year 
is inversely proportional to the magnitude of the average WIMP velocity with respect to 
Earth, — Viab- The largest uncertainty in Viab stems from the Galactic rotation speed at 
the position of the Sun, VcaiRot- Thus measuring 7^ would allow to determine VcaiRot- The 
dependence of 7^ on VcaiRot is shown in Fig. 7. Depending on the value of VcaiRot) between 
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180 km/s and 312 km/s, 7^ could be between 18° and 11°. In Section 4, we estimated that 
between 1400 and a few thousand events would be needed to measure the angular difference 
in recoil direction integrated over six or three months. We used as the basis for our estimate, 
the minimum number of events given in Refs. [15] and [16] needed to determine the mean 
recoil direction with an error of a few degrees. The error in the determination of VcaiRot due 
to an uncertainty in 7^ can be read off of Fig. 7. For example, with an error of 3.5°, the 
measured VcaiRot would be I8OI3Q km/s or 312^79'^ km/s using the three month averaged 
curve, and ISOtlo ^^/^ 3121^90 km/s using the six month averaged curve. 

In Section 5 we showed that with more than a thousand events, it would also be possible 
to detect the annual modulation due to aberration of the energy integrated rate. The inte- 
grated counts over Galactic hemispheres separated by planes perpendicular to Earth's orbit 
would modulate annually, resulting in Galactic Hemisphere Annual Modulations (GHAM), 
with amplitudes larger than the usual non-directional annual modulation. The Galactic 
hemispheres with the largest GHAM are those divided by planes perpendicular to the Earth's 
orbital velocity when T^ab is maximum (see Figs. 8, 9). In this case the full orbital velocity 
is in the direction of one hemisphere at one time and away from it six months later. In the 
examples shown in Figs. 8, 9 and 10, the amplitudes of the GHAM are larger than those of 
the usual annual modulation of the non-directional rate by factor of at least two, and for the 
energy differential rate change sign at much lower energies than the usual non-directional 
annual modulation amplitude. As shown in Fig. 13, assuming no background and perfect 
energy and angular experimental resolution, at least a few ton-yr of S would be required 
to detect the GHAM of the spin-independent energy-integrated rate. To detect the largest 
GHAM of the spin-dependent energy- integrated rate in CF4, a minimum exposure of 60 kg-yr 
would be required as shown in Fig. 14. b. 

The GHAM amplitudes as functions of energy depend strongly on the average velocity 
and velocity dispersion (see Figs. 8, 9 and 10). Work beyond the scope of this paper is nec- 
essary to determine how best aberration features could be used to extract the characteristics 
of the local WIMP velocity distribution. 
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